Theoretical calculations of structural, electronic, and elastic properties of CdSe1−xTex: A first principles study
Shakil M1, †, , Zafar Muhammad2, Ahmed Shabbir2, Hashmi Muhammad Raza-ur-rehman2, Choudhary M A2, Iqbal T1
Department of Physics, Hafiz Hayat Campus, University of Gujrat, Gujrat 50700, Pakistan
Simulation Laboratory, Department of Physics, the Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan

 

† Corresponding author. E-mail: shakil101@yahoo.com; m.shakil@uog.edu.pk

Abstract
Abstract

The plane wave pseudo-potential method was used to investigate the structural, electronic, and elastic properties of CdSe1−xTex in the zinc blende phase. It is observed that the electronic properties are improved considerably by using LDA+U as compared to the LDA approach. The calculated lattice constants and bulk moduli are also comparable to the experimental results. The cohesive energies for pure CdSe and CdTe binary and their mixed alloys are calculated. The second-order elastic constants are also calculated by the Lagrangian theory of elasticity. The elastic properties show that the studied material has a ductile nature.

1. Introduction

One of the basic features of the solar cell is its ability to absorb effectively a wide spectrum of photons contained in the solar radiation. This feature depends on the intrinsic optical and electronic properties of the semiconductor material, and one critical parameter is the energy band gap. The II–VI chalcogenides and their mixed ternary alloys show semiconducting properties.[1,2] These chalcogens and their alloys are especially suitable for the conversion of solar energy to electrical energy in photovoltaic or photo electrochemical devices. For a selected composition range of the solid solution of CdSe and CdTe, the band-gap and lattice parameters can be changed by altering the relative amount of chalcogenides in the CdSexTe1 − x compound. For thin films of CdSe1−xTex, the crystal structure and band-gap can be adjusted by changing the concentrations of Se and Te. The band gap of this material can be varied from 1.39 eV to 1.72 eV by appropriately adjusting the atomic composition (x) from 0 to 1.[3] This range of band gap is important for photovoltaic conversion, because the ideal band gap for the solar cell is 1.4 eV. Therefore CdSe1−xTex ternary thin films are widely used in solar cells,[4,5] solar control, and photovoltaic applications.[6,7] They can also be used for the photo-assisted decomposition of water,[8] etc. Normally, CdSe and CdTe exhibit hexagonal wurtzite and cubic zinc blende (ZB) structures, respectively. A structural phase transition is expected in the CdSe1−xTex pseudo-binary system as a function of the composition. A number of investigations have reported that these films do not have a single phase over the entire composition range and the presence of either the cubic or hexagonal phase depends on the proportion of selenium and tellurium in the pseudo-binary system. Murali and Jayasutha reported the hexagonal phase for the CdSe1−xTex films deposited on a substrate by brush plating at different temperatures.[9] However, a few experimental results were inconsistent with the above discussed results. Uthanna and Reddy,[10] Mangalhara et al.,[11] and Islam et al.[12] have observed the presence of only the cubic phase in the CdSe1−xTex films over the entire composition range.

A large number of theoretical and experimental results[112] have been reported for CdSe1−xTex. However, all the studies focused on pure CdSe and CdTe. In this study, the electronic properties of ternary alloys (i.e. CdSe1−xTex) are calculated by the LDA+U approach based on density functional theory (DFT) and the elastic properties are reported for the first time. It is envisaged that the elastic properties have significant influences on the devices.

2. Computational details

To simulate the various physical properties of CdSe1−xTex alloys, we employed the plane wave pseudo potential method to solve the Kohn–Sham equations. The well known and versatile DFT based Quantum ESPRESSO code[13] was used to study the properties within the framework of DFT[14] and this method produces consistent results for the physical properties of different materials. For the exchange–correlation (XC) potential, the local density approximation (LDA)[15] of Perdew and Zunger (pz)[16] was applied. The ultrasoft pseudo-potential (USPP) was used to describe the interaction between the ionic cores and the valance electrons. Plane waves with a cutoff energy of 30 Ry were used to characterize the electronic wave functions, and the Brillouin zone integrations were performed using an 8 × 8 × 8 Monkhorst–Pack[17] grid of k-point mesh. In the present study, CdSe1−xTex was modeled at various compositions with a step of 0.25 using eight atoms per unit cell. A partial substitution of an element by another in a compound can cause drastic changes in the physical properties of the parent material. In this study, CdSe is the parent material and Se is partially substituted by Te. The total energy for each concentration was calculated and it was found that CdSe and CdTe are more stable in the zinc-blende structure. The equilibrium lattice parameters and bulk modulli of CdSe1−xTex were obtained by calculating the total energy at different lattice constants and fitting the data to Murnaghan’s equation of states. As LDA and GGA both calculated a very small energy gap, a Hubbard U correction[18] was incorporated with LDA for the determination of the electronic band structure of CdSe1−xTex. The elastic constants and moduli were calculated by the Lagrangian theory of elasticity.[19]

3. Results and discussion
3.1. Structural properties

As a starting step, the total energy was calculated as a function of the lattice constant for binary compounds CdSe, CdTe and their mixed alloys. The energy versus lattice constant data are plotted in Fig. 1. The data were fitted to Murnaghan’s equation of state[20] to determine the equilibrium lattice constant a0 and the bulk modulus B. The optimized results are summarized in Table 1 to compare with other theoretical and experimental data. The optimized results for a0 and B are presented in Fig. 2 along with Vegard’s law. To study the alloys, it is assumed that all atoms in the crystals are located at their ideal sites. It is also assumed, according to Vegard’s law, that any compound has a linear variation in lattice constant with composition x.[21] Our calculated lattice parameters at different compositions for CdSe1−xTex alloys are found to vary exactly linearly, consisting with the law. The calculated lattice parameters of the CdX compound increase in the following sequence: a0(CdSe) < a0(CdSe0.75Te0.25) < a0(CdSe0.50Te0.50) < a0(CdSe0.25Te0.75) < a0(CdTe). The reason is that the tetrahedral covalent radius of Te (0.134 nm) is larger than that of Se (0.115 nm).[22] The cohesive energy is a measure of the strength of the forces which for the ZB structure of CdX is calculated as

where is the total energy of CdSe and CdTe with equilibrium lattice constants, and and are the energies of the isolated Cd and Se/Te atoms in the free state, respectively. The calculated cohesive energies, listed in Table 1, are very close to the experimental values.[23] The accuracy of the cohesive energies calculated by Quantum ESPRESSO is encouraging as shown by this comparatively satisfactory result, confirming that the calculated results in this work are reliable. Further analysis of these results shows that CdSe has relatively higher structural stability than CdTe due to the larger absolute value of the cohesive energy.

Fig. 1. Structural optimization plots for CdSe1−xTex alloys in the ZB phase: (a) x = 0.00, (b) x = 0.25, (c) x = 0.50, (d) x = 0.75, (e) x = 1.00.
Fig. 2. Variations of the calculated (a) lattice constants and (b) bulk moduli with x for CdSe1−xTex alloys.
Table 1.

Calculated lattice parameters and bulk moduli of CdSe1−xTex alloys compared with experimental and other theoretical results.

.
3.2. Electronic properties

The applications of the studied materials require precise knowledge of the energy band gaps along with the aligned placement of main conduction band (CB) valleys. Band gaps and band structures provide esteemed information for the fabrication and development of electronic, opto-electronic, and magneto-electronic devices. To explore the electronic properties of CdSe, CdTe compounds and their alloys, we calculated the band structure at high symmetry points in the first Brillouin zone in the ZB structure within LDA and LDA+U at equilibrium lattice constants. As the valence band maximum (VBM) and the conduction band minimum (CBM) locate at the same point, CdSe1−xTex has a direct energy band gap. This direct energy band gap decreases from CdSe (x = 0) to CdTe (x = 1). The obtained results of this study are compared with the experimental and other theoretical results in Table 2. In the case of LDA, the band gaps are smaller than the experimental values, since the LDA approach underestimates the band gaps for all materials of interest. While in the case of LDA+U, the energy band gaps are very close to the experimental results. The value of U (in eV) decreases from CdSe (x = 0) to CdTe (x = 1) as listed in Table 2. The band structures obtained with LDA and LDA+U are depicted in Fig. 3. In this figure, we notice that the band gap varies linearly with the composition (0.25, 0.5, and 0.75). Figure 3 also shows that the bottom of the conduction bands and the top of the valence bands are situated at the Γ point in the Brillouin zone, resulting in a direct band gap. To elucidate the nature of the electronic band structure, we calculated the total and partial density of states (PDOS) of the CdSe1−xTex alloys and present the results in Fig. 4. The first part VB1 is due to the s states of Se and Te. The second part VB2 is due to the 3d states of Cd. Whereas the uppermost valence band, VB3, is formed due to Se and Te p states along with Cd s states between −6 eV and 0 eV except for CdSe where this range is from −4 eV to 0 eV. While the Se/Te s states (VB1) are in the range from −12.0 eV to −10.0 eV. Similarly, the lowest CB is dominated by Cd s and p chalcogen states. It can be seen from the plot that the Te s state is closer to the Cd d state as compared to the Se s state. It can be observed clearly from Fig. 4 that when x changes from 0.0 to 1.00, the density of states (DOS) is changed. In general, the band gap of CdTe is smaller than that of CdSe, but in our case (by LDA), it is larger as shown in Table 2. Similar behavior has been reported.[28,29] It is recovered with the help of the LDA+U approach. It confirms that the Te incorporation into CdSe induces an obvious band gap narrowing. It is necessary to point out that the VBM is always in the Fermi level and it has no obvious shift, which can be clearly seen in Fig. 3.

Fig. 3. Electronic band structures of CdSe1−xTex: (a) x = 0.00, (b) x = 0.25, (c) x = 0.50, (d) x = 0.75, (e) x = 1.00.
Fig. 4. Total and partial DOSs of CdSe1−xTex: (a) x = 0.00, (b) x = 0.25, (c) x = 0.50, (d) x = 0.75, (e) x = 1.00.
Table 2.

Comparison of calculated band gaps of CdSe1−xTex alloys with experimental and theoretical results in the ZB structure.

.
Fig. 5. Variation of energy gap of CdSe1−xTex with the concentration of Te.
3.3. Elastic properties

The elastic constants are important and necessary parameters to describe the response to an applied macroscopic stress and especially important because they are correlated to numerous physical properties, such as the elasticity, mechanical stability, and stiffness of the materials. The cubic crystals have only three independent elastic constants, i.e., C11, C12, and C44. They can be determined by computing the stress generated by applying a small strain to an optimized unit cell. For three structures, we carried out calculations with strains in the range of −0.05 to 0.05 for each distortion. The systems were fully relaxed after each distortion in order to reach the equilibrium state with approximately zero forces on all atoms. In this work, we calculated the second-order elastic constants by the analysis of change in calculated stress resulting from change in the strain. The obtained results are given in Table 3, which show good agreement with other theoretical and available experimental results. The deliberated results about bulk moduli by energy minimization and from the elastic constants are comparable. It confirms the accuracy and consistency of our calculated elastic constants. The bulk modulus of CdTe is smaller than that of CdSe, which approves the legitimacy of the known empirical relation proposed by Cohen[33] B = 1761 × d−3.5, where d is the nearest-neighbor distance in angstroms. The traditional and well known mechanical stability conditions of elastic constants for cubic crystals are C11C12 > 0, C11 > 0, C44 > 0, C11 + 2C12 > 0, and C12 < B < C11.[34] For CdSe1−xTex semiconductor alloys, the obtained elastic constants are given in Table 3 and plotted in Fig. 6 which satisfy the above stability criteria, and thus corroborate the mechanical stability of the structure. Since the elastic anisotropy of materials has an imperative application in engineering, it is highly associated with the possibility to induce micro-cracks in the material. To enumerate the elastic anisotropy of the alloy, we computed the anisotropy factor A from the elastic constants using the following expression:

For an absolutely isotropic material, A is equal to 1, while any value smaller or larger than 1 indicates an anisotropic material. The magnitude of deviancy from 1 is a measure of the degree of elastic anisotropy possessed by the crystal. From the computed anisotropy factors given in the Table 3, one can see that A is greater than 1.0 and we cannot consider CdSe1−xTex as elastically isotropic. Once the elastic constants C11, C12, and C44 are calculated, the related properties, namely, the shear modulus G, the Young modulus E, and the Poisson ratio v can be estimated under Voigt–Reuss–Hill approximations[35] using the following relations:

where B = (C11 + 2C12)/3 is the bulk modulus, and G = (GR + GV)/2 is the isotropic shear modulus. Here GR is the Reuss shear modulus corresponding to the lower bound of G, and GV is the Voigt shear-modulus corresponding to the upper bound of G; they can be written as GV = (C11C12 + 3C44)/5 and 5/GR = 4/(C11C12) + 3/C44. The calculated elastic moduli for polycrystalline CdSexTe1 − x are also given in Table 3.

Table 3.

Comparison of calculated elastic constants of CdSe1−xTex alloys with other results in ZB structure.

.

The shear modulus is related to the resistance to plastic deformation and the bulk modulus is associated with the resistance to fracture. Pugh[36] introduced the new criteria about the brittleness and hardness of a material which is related to the ratio of the bulk modulus to the shear modulus (B/G). A high (low) value of B/G means a ductile (brittle) material. The critical value which separates ductility from brittleness is about 1.75. The results infer that none of the CdSe1−xTex compounds are brittle. Poisson’s ratio v provides further information for dealing with the characteristic of the bonding forces. The values 0.5 and 0.25 are the upper and lower boundaries for central force solids, respectively.[37] Since the measured Poisson ratio is between these values, so the inter-atomic forces are central in these compounds.

Fig. 6. Variation of (a) elastic constants and (b) elastic modulus in CdSe1−xTex at various concentrations of tellurium in ordered alloys.
4. Conclusions

Employing the pseudo-potential technique based on DFT within LDA and LDA+U, we have investigated the structural, electronic, and elastic properties of the binary and ternary alloy of CdSe and CdTe. The material was chosen because of its vast field of application. This study reveals the composition dependent changes in the ground state properties such as lattice constants (a0), bulk moduli, and elastic properties of CdSe1−xTex semiconductor alloy. The calculated results are in good agreement with the theoretical and experimental results. This study gives confirmations of the LDA+U approach that was used to improve the underestimated electronic band gaps of Cd-based selenium, tellurium and ternary semiconductor alloys. The simulated band structures and DOSs show that the CdSe1−xTex alloys are semiconductor in nature with a small direct band gap. For this calculation, strong intra-correlation factor U for Cd 4d electrons was used and the results with the LDA+U approach are comparable to the experimental results. No deviation in lattice constant from Vegard’s law is observed for CdSe1−xTex. A linear dependence on composition x is observed for the bulk modulus and lattice constant of the studied alloys. The B/G value reflects that the studied material has a ductile nature for the whole concentration.

Reference
1Zafar MAhmed SShakil MChoudhary M A 2014 Chin. Phys. 23 106108
2Zafar MAhmed SShakil MChoudhary M AMahmood K 2015 Chin. Phys. 24 076106
3Muthukumarasamy NBalasundaraprabhu RJayakumar S 2004 Physica Status Solidi a 201 2312
4Hodes GManassen JCahen D 1980 J. Am. Chem. Soc. 102 5962
5Muthukumarasamy NVelumani SBalasundaraprabhu RJayakumar SKannan M D 2010 Vacuum 84 1216
6Sebastian P JSivaramakrishnan V 1990 J. Phys. D: Appl. Phys. 23 1114
7Kumar LSingh B PMisra AMisra S C K 2005 Physica 363 102
8Velumani SMathew XSebastian P J 2003 Solar Energy Materials and Solar Cells 76 359
9Murali K RJayasuthaa B 2009 Solar Energy 83 891
10Uthanna SReddy P J 1983 Solid State Commun. 45 979
11Mangalhara J PThangaraj RAgnihotri O P 1989 Solar Energy Materials 19 157
12Islam RBanerjee H DRao D R 1995 Thin Solid Films 266 215
13Paolo GStefano B 2009 J. Phys.: Condens. Matter 21 395502
14Hohenberg PKohn W 1964 Phys. Rev. 136 864
15Perdew J PWang Y 1992 Phys. Rev. 45 13244
16Perdew J PZunger A 1981 Phys. Rev. 23 5048
17Monkhorst H JPack J D 1976 Phys. Rev. 13 5188
18Anisimov V IZaanen JAndersen O K 1991 Phys. Rev. 44 943
19Wallace D C1998Thermodynamics of Crystals, Dover Books on PhysicsNew YorkDover Publications
20Murnaghan F D 1944 Proc. Natl. Acad. Sci. USA 30 44
21Vegard L 1921 Zeitschrift für Physik 5 17
22Pyykkö P 2012 Phys. Rev. 85 024115
23Madelung OSchultz MWeiss H1982Landolt-Bornstein Numerical Data and Functional Relationships in Science and TechnologyBerlinSpringer-Verlag
24Ameri MFodil M 2012 Materials Sciences & Applications 3 768
25Cardona MMadelung OHabeke GRössler U1982Semiconductors Physics of Group IV Elements and III–V CompoundsHellwege K HMadelung OLandolt-Bornstein, New series, Group III17
26Al-Douri YReshak A HBaaziz HCharifi ZKhenata RAhmad SHashim U 2010 Solar Energy 84 1979
27Guo LZhang SFeng WHu GLi W 2013 J. Alloys Compd. 579 583
28Ouendadji SGhemid SMeradji HHassan F E H 2010 Computational Materials Science 48 206
29Fleszar AHanke W 2005 Phys. Rev. 71 045207
30Khan IAhmad IAliabad HMaqbool MComputational Physics(physicscompph) arXiv:120108702012
31Rozale HLazreg AChahed ARuterana P 2009 Superlattices and Microstructures 46 554
32Patil SPawar P2012Biological and Physical Sciences21472
33Cohen M L 1985 Phys. Rev. 32 7988
34Shein I RIvanovskii A L 2008 J. Physics: Condens. Matter 20 415218
35Chung DBuessem W 1967 J. Appl. Phys. 38 2535
36Pugh S F 1954 The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 45 823
37Dao MChollacoop NVan Vliet KVenkatesh TSuresh S 2001 Acta Materialia 49 3899
38Sharma SVerma A SSarkar B KBhandari RJindal V K2011AIP Conference Proceedings139322910.1063/1.3653693
39Deligoz EColakoglu KCiftci Y 2006 Physica 373 124
40Kheloufi NBouzid A 2016 J. Alloys Compd. 659 295